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ABSTRACT 

The dynamical mass of clusters of galaxies, calculated in terms of modified Newto- 
nian dynamics, is a factor of two or three times smaller than the Newtonian dynamical 
mass but remains significantly larger than the observed baryonic mass in the form of 
hot gas and stars in galaxies. Here I consider further the suggestion that the unde- 
tected matter might be in the form of cosmological neutrinos with mass on the order 
of 2 eV. If the neutrinos and baryons have comparable velocity dispersions and if the 
two components maintain their cosmological density ratio, then the electron density 
in the cores of clusters should be proportional to T 3 / 2 , as appears to be true in non- 
cooling flow clusters. This is equivalent to the "entropy floor" proposed to explain the 
steepness of the observed luminosity-temperature relation, but here preheating of the 
medium is not required. Two fluid (neutrino-baryon) hydrostatic models of clusters, 
in the context of MOND, reproduce the observed luminosity-temperature relation of 
clusters. If the (3 law is imposed the gas density distribution, then the self-consistent 
models predict the general form of the observed temperature profile in both cooling 
and non-cooling flow clusters. 



1 INTRODUCTION 



In the past 25 years, Milgrom's proposed alternative to 
dark matter, the modified Newtonian dynamics or MOND 
(Milgrom 1983), has enjoyed considerable phenomeno- 
logical success with respect to galaxy scaling relations 
and individual rotation curves (see Sanders & McGaugh 
2002 for a review). The range of this predictive power 
has recently been extended from very low luminosity 
objects (Milgrom & Sanders 2007 1 to luminous early-type 
spirals (Sande rs fc Noordermeer 2007[) . However, for some 
time it has been recognised that MOND does not fully 
account for the mass discrepancy in rich clusters of galaxies 
l|The fc White 19881 IGerbal et al. 19921 ISanders 19991 
Aguirre, Schaye fc Quataert 2001] ISanders 20031 

Pontecouteau fc Silk 2005^ This is basically because, 
with MOND, the departure from Newtonian dynamics 
appears below a critical acceleration , ao- The value of 
oo determined from galaxy rotation curves is about 10 -8 
cm/s 2 . In clusters of galaxies, the observed acceleration 
(estimated from the hydrostatic gas equation) is typically 
greater than ao in the central regions, so MOND cannot 
alleviate the observed discrepancy there. When the MOND 
version of the hydrostatic gas equation is applied to the 
observed temperature and density distributions in X-ray 
emitting clusters of galaxies, one finds that there remains a 
discrepancy roughly of a factor of three between the iden- 
tifiable baryonic mass in gas and stars and the dynamical 
mass. That is to say, MOND also requires the presence of 
undetected, or dark, matter in clusters of galaxies, albeit 
significantly less than conventional Newtonian dynamics. 
X-ray and weak lensing observations of the famous "bul- 



let cluster" (fCTowc ct al. 2006} are now presented as defini- 
tive evidence for non-dissipative dark matter in clusters of 
galaxies and, by extension, as evidence against MOND. Had 
it been previously claimed that MOND fully accounts for the 
kinematic and X-ray observations of clusters without dark 
matter, then this object would indeed have been quite prob- 
lematic for the theory. But, in fact, the bullet creates no 
additional difficulties for MOND; the quantity of dark mat- 
ter required is consistent with that suggested by the previ- 
ous analyses (Angus et al. 2007). What the bullet cluster 
adds (and this is a significant addition) is convincing evi- 
dence that the dark component cannot be dissipative like 
the extended X-ray emitting gas. 

At first thought, it might appear negative for MOND 
that additional dark matter is required in the cluster en- 
vironment. But one could also consider this to be a bold 
prediction- that there is more matter to be discovered in 
clusters. For example, there are more than enough unde- 
tected baryons to make up the missing dark component; they 
need only be present in some non-dissipative form which 
is difficult to observe. Moreover, it is now known that at 
least two of the three active neutrino types have non-zero 
mass (Fukuda ct al. 1998). Primordial neutrinos are present 
in the Universe with a number density comparable to that 
of photons, so non-baryonic dark matter certainly exists. If 
the neutrino mass scale is as large as 2 eV, then neutri- 
nos would comprise a non-baryonic component of rich clus- 
ters (Sanders 2003), but because of phase space constraints 
l|Trem ainc & Gunn 1979) they could not accumulate in in- 
dividual galaxies. (3 decay experiments have restricted the 
mass of the electron neutrino to be less than about 2.2 eV, 
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so it remains possible that this is the appropriate mass scale 
l|Groom et al. 2000)1 . The addition of cosmological neutri- 
nos with mass of this order produces a reasonable fit to the 
observed angular power spectrum of the cosmic microwave 
background radiation, through the second peak, without the 
presence of cold dark matter (McGaugh 2004} . 

Therefore, here I consider further the possibility that 
the missing component in clusters of galaxies in the con- 
text of MOND may be neutrinos. I demonstrate that the 
form of correlations between observable quantities- the size- 
temperature relation (R oc T) and the gas mass-temperature 
relation (M g oc T 2 ) are consistent with MOND expectations 
without any additional astrophysical mechanism other than 
violent relaxation. However, the luminosity-temperature re- 
lation (L oc T 3 ) is not consistent with MOND by it- 
self; an additional dark component seems to be required 
IjSanders 2003[> . 

For a sample of well-studied clusters, the MOND dy- 
namical mass is proportional to but a factor of three or four 
times larger than the observable gas mass. This suggests 
that the detectable baryonic mass is a fixed fraction of the 
dark mass. If the mass of each active neutrino type is about 
2 eV, then the ratio of neutrino to baryonic cosmological 
mass densities would be about 2.8, consistent with the in- 
ferred (via MOND) ratio of gas to dark matter in clusters. 

I consider the structure of two-component isothermal 
objects consisting of neutrinos and gas. The neutrino fluid 
is described by the equation of state of partially degenerate 
fermions, where the degree of degeneracy in the centre is 
arbitrarily set to a specified value. Because the maximum 
density of neutrinos is proportional to J 13 / 2 this means that, 
in a uniformly mixed two-component fluid, the central gas 
density would also be proportional to T 3//2 . This appears 
to be the case in non-cooling flow clusters where there has 
been no inward flow of gas resulting in a rearrangement 
of the cosmological density ratio of baryons-to-neutrinos. 
With such models, the observed luminosity- temperature 
scaling relation for clusters is recovered primarily because 
L oc n 2 oc T 3 . 

I then consider individual clusters with the two-fluid 
model. After specifying the mass scale of the neutrinos and 
the degree of degeneracy in the centre, I assume that the 
presumably constant neutrino velocity dispersion is equal 
to that of the gas as implied by the mean emission-weighted 
gas temperature. Given the observed gas density distribu- 
tion (via /3-model fits to the X-ray intensity distribution) 
I calculate the radial dependence of the gas temperature 
which is consistent with this density distribution. The only 
free parameter is the central gas temperature which must lie 
in a narrow range in order to yield a sensible temperature 
distribution (one in which the temperature does not rapidly 
increase to large values or fall to zero) . These calculated tem- 
perature distributions are found to be generally consistent 
with those now observed- a slowly decreasing temperature 
for non-cooling flow clusters, and a temperature which first 
increases with radius and then decreases for cooling flow 
clusters. 

I conclude that the two-fluid neutrino-baryon model for 
clusters of galaxies in the context of MOND, is consistent 
both with the correlations between observable quantities and 
with the observed density and temperature distributions in 
individual clusters. 



In all that follows I have scaled observational results 
and correlations to Ho = 72 km/s-Mpc. 



2 GLOBAL SCALING RELATIONS FOR 
CLUSTERS OF GALAXIES 

It is of interest to consider correlations between the ob- 
served, and not the inferred, properties of X-ray emitting 
clusters of galaxies. One of the most obvious correlations is 
that between the cluster size (out to a specified X-ray inten- 
sity level) and the temperature of the hot gas. Mohr et al. 
(2000) find that this relation, for nearby clusters, is: 



R = 0.5 



T 



6keV 



Mpc 



(1) 



with a surprisingly low scatter of about 15%. The expec- 
tation for self-similar collapse is more like R oc T a where 
a ~ 1/2 — 2/3. Mohr et al. explain the difference between 
the observations and the expectations by noting that the 
fraction of gas (as opposed to the total baryonic fraction 
including the luminous matter in galaxies) appears to in- 
crease with cluster temperature. While this seems to be a 
general trend (|David et al. 1990l[Edge fc Stewart 1991} , the 
observed relationship between gas mass/total baryonic mass 
and temperature exhibits enormous scatter. It is unclear how 
such a tight correlation can survive. 

In terms of MOND, eq. 1 takes on a different mean- 
ing. The mean internal gravitational acceleration in clusters 
is given roughly by a = a 2 /R (a is the velocity disper- 
sion and proportional to temperature). Then we may in- 
terpret eq. 1 as defining a single characteristic internal ac- 
celeration which is independent of cluster size or temper- 
ature; i.e., a oc T/R « 6 x 10 -9 cm/s 2 . This is within 
a factor of two of the MOND critical acceleration (do ~ 
10~ 8 cm/s 2 ). It has been pointed out previously that in 
pressure supported, nearly isothermal systems, the inter- 
nal acceleration is approximately the MOND acceleration 
(Sanders & McGaugh 2002). Viewed in this way, the clus- 
ter temperature-size relation simply reflects that character- 
istic internal acceleration and requires no astrophysical in- 
put other than the requirement of a near isothermal state, 
presumably due to violent relaxation. 

A second correlation is that between the temperature 
and the observed gas mass. Mohr et al. (1999) give 
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For self-similar collapse, making use of the Newtonian virial 
theorem, the expectation would be M g oc T 3 / 2 . The dis- 
agreement between observations and expectations is again 
attributed to a systematically changing fraction of baryons 
in hot gas perhaps resulting from increasing energy injec- 
tion into the intra-cluster medium or a decreased efficiency 
of galaxy formation in hotter clusters. 

With MOND, however, eq. 2 would be the expecta- 
tion from the dynamics of isothermal (or semi-isothermal) 
spheres; it is, in effect, an extrapolation of the Faber- Jackson 
relation for elliptical galaxies ([Sanders 1994]) and, in the 
case of MOND would apply to all pressure supported, near 
isothermal systems ( |Milgrom 1984} . Basically, the equation 



Neutrinos 3 



of hydrostatic equilibrium in the MOND regime yields 
M 



•06 T kev 



(3) 



1O 14 M 

for the total dynamical mass in terms of MOND. Here I have 
assumed an isothermal /3-model with an outer logarithmic 
density gradient of —3/3 where j3 = 0.6 on average. Although 
the form of the temperature-gas mass relation is predicted, 
the normalisation is not; the MOND dynamical mass implied 
by eq. 3 is roughly 3.6 times larger than the observed gas 
mass. 

The third correlation, and the most discussed in the lit- 
erature, is that between luminosity and temperature. Again 
the correlation is well-fitted by a power-law; for example, 
Ikebe et al. (2002) give 



■■ 4 x 10 42 T fceV - 2 - 5 ergs/s 2 



(4) 



for the total X-ray luminosity between energies of 0.1 to 
2.4 keV. Others find exponents closer to 3 with suggestion 
of a steepening for cooler clusters (|Arna ud & Ev rard 1999[) . 
The expectation for self similar cluster formation in the con- 
text of CDM is a shallower power law L oc T 2 . Here, pre- 
heating of the intergalactic medium (heating before clus- 
ter formation) by winds from forming galaxies is invoked 
to to provide an "entropy floor" for the subsequent intra- 
cluster medium (entropy is defined as Tk e v /ra e 2//3 ). This 
would have the effect of inflating the the intra-cluster 
medium particularly in lower temperature galaxies and 
steepening the predicted luminosity-temperature relation 
(Ponmon, Cannon & Navarro 1999). 

MOND, in this case, fairs worse. In MOND a virial- 
ized system has a characteristic density pM oc ao 2 /(GT). 
As we saw above, the characteristic radius of the cluster is 
R oc T/ao- For free-free radiation, L oc p T 1/2 R 3 ; therefore, 
the naive prediction would be L oc T 3 / 2 . I discussed this 
point previously (Sanders 2003) and noted that a dark mat- 
ter component with a constant density and a core radius 
roughly twice that of the gas core radius could bring the 
MOND prediction in line with the observations. 



3 THE MAGNITUDE OF THE DISCREPANCY 
WITH MOND 

MOND predicts more mass than is directly observed in clus- 
ters of galaxies. Determination of the magnitude of this dis- 
crepancy, however, is not straightforward. For an isothermal 
gas, the hydrostatic gas equation in the MOND limit implies 
that the mass within radius R is given by 

2 (^) 2 (Ga )- 1 (5) 
\wm p J 



where w = 0.62 is the mean atomic weight, m p is the pro- 
ton mass, and a = din pj din r is the logarithmic density 
gradient evaluated at R. The gas density distribution is typ- 
ically described by the /3 model used to fit the X-ray surface 
brightness distribution: 
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Therefore beyond three or four core radii a — > —3/3- i.e., the 
predicted dynamical mass converges. On the other hand, 
= 0.5 — 1.0 which means that the gas mass continues 
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Figure 1. The logarithm of the MOND dynamical mass plotted 
against the logarithm of inferred baryonic mass (both in units of 
10 14 Mq). The crosses indicate cooling flow clusters (those where 
the central cooling timescale is less than the Hubble time), and 
the solid points are non-cooling flow clusters. 



to increase with radius (the /3 model surely must have a 
limited range of viability). Therefore, the magnitude of the 
discrepancy (defined as Md/M g ), as estimated from the /3 
model, decreases with radius; i.e., the discrepancy depends 
upon the radius within which the estimate is made. 

I have determined the MOND dynamical mass and gas 
mass for a well-studied sample of clusters given by Reiprech 
& Bohringer (2002). I make use of the radius-temperature 
relation to estimate the cutoff radius, but I take this limiting 
radius to be about 20% larger than that given by eq. 1. This 
implies that the cutoff is, on average, about six core radii 
which is a probable range of validity for the /3 model fit. 
Since the accelerations in the central regions of the clusters 
are typically on the order of ao, the form of the MOND 
interpolating function plays a role; here I take the simple 
form suggested by Zhao and Famaey (2005): 



p(x) = 



(7) 



1 +x 

where x = a/a a with a = akT/(wm p R). The MOND dy- 
namical mass is then given by 

xM N 



M„ 



1 +X 



(8) 



where the Newtonian mass is Mn = —aRkT / (Gwm p ) . The 
gas is only a fraction of the total baryonic mass- the re- 
mainder being in the form of stars in galaxies. As before, 
I take the ratio of intra-cluster gas to stellar mass as given 
by Tkev /2.5 ([Sanders 1999). There is considerable scatter 
about this approximate relation. 

The result of this analysis is shown in Fig. 1 which is 
a plot of the MOND dynamical mass against the inferred 
baryonic mass out to a limiting radius given by 100Ti; e y 
kpc. The ratio of the MOND dynamical mass to the bary- 
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onic mass is, on average, 3.7 with no significant difference 
between the cooiing flow and non-cooling flow galaxies. That 
is to say, the observed baryonic mass in clusters accounts for 
one-third to one-forth of the MOND dynamical mass. 

It is evident from Fig. 1 that the ratio of dynamical to 
baryonic mass is roughly independent of the baryonic mass. 
There is an indication of a deviation at low values of the 
baryonic mass possibly because of an increased fraction of 
matter in the stellar component of galaxies, but overall, the 
points lie parallel to the line of equality. This is of inter- 
est because it suggests that a fixed cosmic ratio of dark to 
baryonic matter may be sampled by the cluster medium. 



4 THE ROLE OF MASSIVE COSMOLOGICAL 
NEUTRINOS 

If neutrinos are more massive than a few tenths of an elec- 
tron volt, then, because of the small mass differences, all 
types have about the same mass. The cosmological density 
of neutrinos would then be given by f2„ = 0.062m„ where 
m v is the mass of a single neutrino type in eV. In this re- 
lation I have assumed three left-handed neutrino types and 
their anti-particles, but the formula would apply for Majo- 
rana as well as Dirac neutrinos. 

Given that the cosmological density of baryons is fij, = 
0.044 (e.g. Spergel et al. 2006, with h=0.72), this means that 

flv/Qi = 1.4m„. (9) 

If we suppose that clusters of galaxies sample this universal 
neutrino- baryon density ratio then the dynamical mass-to- 
baryonic mass ratio implied by Fig. 1 would mean that m v w 
2 eV. 

Cosmological neutrinos freeze out of the primordial fire- 
ball below temperatures characteristic of the weak interac- 
tion scale of a few meV. The phase space distribution is 
Fermi-Dirac with a maximum for each species of one parti- 
cle per cell with volume 2h 3 ; i.e., one-half that of complete 
degeneracy. This initial phase space density is maintained 
as a limit on the final phase space density of any collapsed 
virialized object (Tremaine & Gunn 1979); of course, the 
degeneracy limit is absolute. 

In the formation of a cluster scale object out of the 
mixed fluid of neutrinos and baryons, it is expected that the 
two fluids attain the same velocity dispersion via violent re- 
laxation (Kull, Treumann & Bohringer 1997). With this as- 
sumption, the maximum contribution of each neutrino type 
(including its anti-particle) to the density of the neutrino 
fluid is 



Pu = (2tt) 
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( |Kull, Treumann fc Bohringcr ~1996[ ). For three neutrino 
types (and their anti-particles) this translates to a mass den- 
sity of 
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If the ratio of the densities of baryons and neutrinos in clus- 
ters is the same as the cosmological ratio (eq. 9), this means 
that the density of electrons, in a fully ionised plasma, would 
be 
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Figure 2. The logarithm of of the central electron density (from 
the fitted (3 models of Reiprech & Bohringer) vs. the logarithm 
of the temperature in eV. The points are the non-cooling flow 
clusters and the crosses are the cooling flow clusters. The line 
corresponds to eq. 12 for a neutrino mass of 1.9 eV. 



Therefore, if the ratio of baryon to neutrino densities 
in clusters is identical to the cosmic ratio, as one might ex- 
pect if there has been no subsequent cooling and inflow of 
baryons, then the central electron density should increase 
as T 3 / 2 . This is an observational prediction which can be 
tested, and in Fig. 2 we see the central electron density of in 
the cluster sample of Reiprech & Bohringer plotted against 
the mean electron temperature. The solid points indicate 
those clusters where the cooling time is more than 10 10 years 
(non-cooling flow clusters) and the crosses are those clusters 
with cooling timescales less than this value (cooling flow 
clusters). The solid line is eq. 12 with m v = 1.9 eV. We see 
that in those objects where no cooling and inflow of gas is 
expected- where the cosmic ratio of baryons to neutrinos is 
maintained- there does appear to be such a correlation. 

It is interesting to note that the electron density- 
temperature relation imposed by the neutrino density limit 
and the assumption of fixed baryon-to-neutrino ratio (eq. 12) 
would also correspond to a constant entropy, Tk e v /w e 2//3 . 
Imposing such a temperature dependent limit on the elec- 
tron density is equivalent to imposing an entropy floor. This 
is, in fact, the very mechanism which has been suggested 
as to steepen the predicted luminosity-temperature relation, 
but here the limit is not generated by preheating but by the 
phase space constraints on the neutrino fluid. It is interest- 
ing that the entropy limit corresponding to eq. 12 would 
be 230 keV cm 2 for m„ — 2 eV which is very near the 
value proposed to solve problem of cluster scaling relations 
( Lloyd-Davies, Ponman & Cannon 2000]). 

Before considering the implications for the luminosity- 
temperature relation it is necessary to describe, in a general 
way, the structure of a self-gravitating object consisting en- 
tirely of neutrinos near the degeneracy limit. The object 
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can be approximated as a 7 = 5/3 polytrope; the structure 
would be roughly that of a constant density core with a rapid 
decline beyond a core radius. The core radius can be esti- 
mated by combining the density temperature relation (eq. 
10 or 11) with an appropriate temperature- mass (or viral 
relation). There are two such possible relations. The first, 
for Newtonian dynamics, is 



(13) 



where a is the one-dimensional velocity dispersion, and I 
have taken the potential energy of a 7 = 5/3 polytropic 
sphere. The second is for modified Newtonian dynamics, 
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QMilgrom 1998) . Then with eq. 11 we find two expressions 
for the core radius: Newtonian, 



and MOND, 
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These two values are equal when 
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This meaning of the final expression is that when the tem- 
perature is less than about 3 keV (for m„ w 2 eV) the cluster 
core is in the MOND regime and eq. 16 applies; for higher 
temperatures the cluster core is the Newtonian regime with 
core radius given by eq. 15. In either case, this characteristic 
radius is on the order of 800 kpc and depends very weakly 
on temperature in both regimes. 

Now assuming that baryons are mixed into this cluster 
with the cosmic ratio, most of the X-ray emission would 
be coming from the core region. The X-ray luminosity is 
L oc n e 2 T 1/2 R c 3 or 

LocT 4 

for cooler clusters (MOND regime) and 

l oc r 3 

for warmer clusters (Newtonian regime). In this way, the 
steeper than expected luminosity-temperature relation can 
be understood in terms of the mixed neutrino-baryon fluid. 
The transition from Newton to MOND would explain the 
steepening of the correlation for lower temperature clusters. 

This, of course, is only approximate. X-ray emission also 
originates beyond the core, and the gas is not isothermal. In 
the next section I derive detailed two fluid models of clusters. 



5 MOND-NEUTRINO-BARYON MODELS OF 
CLUSTERS 

The contribution of neutrinos to the mass bud- 
get of clusters has been considered previously 
in the context of mixed CDM-HDM models 
(|Kofman et al. 19961 |Treumann, Kull fc Bohringer 2000| 
|Nakajima fc Morikawa 2 007[ ). In the process of collapse 



and virialization of a multi-fluid mixture of neutri- 
nos and baryons (and/or CDM), violent relaxation 
( |Lynden-Bell 1967| is expected to produce equal velocity 
dispersions for the various components, i.e., Gb = a v . I 
assume this condition here, although it may not be strictly 
true; for example, violent relaxation may be incomplete 
in the outer regions of clusters. Moreover, recent observa- 
tions have demonstrated that the hot gas in clusters is not 
isothermal, in general. The implications for the hypothetical 
neutrino fluid are unclear, but I will assume an isothermal 
condition with the neutrino velocity dispersion equal the 
emission-weighted velocity dispersion of the gas. 

It should also be kept in mind that the density given 
by eq. 11 is an upper limit to the density of the neutrino 
fluid; this value follows from the limit imposed by initial 
phase space density via the collision-less Boltzmann equa- 
tion (Trcmain e fc Gunn 19T9 j) : the course grained phase 
space density in a final relaxed virialized system, and hence 
the space density, may be less significantly less than this 
value. Moreover, while the neutrino fluid may be partially 
degenerate in the core, it will certainly depart from degener- 
acy in the outer regions where the space density declines and 
the velocity dispersion remains constant. Therefore, with re- 
spect to the neutrino fluid, we need to consider the equation 
of state of partially degenerate matter. 

With these caveats in mind, I calculate the structure 
of the neutrino-baryon fluid by applying the hydrostatic gas 
equation separately to each fluid: 



1 dP r 
Pi dr 



-9 



(18) 



where the subscript i refers to one of the two fluid compo- 
nents, Pi is the pressure of that component, pi is the density, 
and g is the total gravitational acceleration resulting from 
the two components and given by the simple MOND expres- 
sion, 



gfi(g/a ) = g N 



(19) 



with p, given by eq. 7 and the gN, the Newtonian force, given 
by g N = G{M v {r) + M b {r))/r 2 

The pressure of the baryonic component is 



D 2 
Pb = Pb&b 



(20) 



The velocity dispersion for the baryons a b 2 (the tempera- 
ture) may be specified as a constant (isothermal) or may 
vary with radius if the density distribution is specified. 

The equation of state for a partially degenerate neutrino 
gas is given parametrically by 



and 
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Pv = o F S(X) 



(21) 



(22) 



where g is the statistical weight, 77 = m^/h 3 , and 

p CO 

F P (X)= / x^l + expix-xT'dx (23) 



(see e.g. Landau & Lifshitz 1980). Here the degeneracy fac- 
tor x is the chemical potential divided by the temperature; 
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log(radius) kpc log T k 



Figure 3. The logarithm of of the neutrino fluid density with 
m„ = 1.9 eV (solid curve) and the baryon fluid density (dashed 
curve) as a function of the logarithm of radius. The two fluids are 
assumed to be isothermal corresponding to the gas temperature of 
5 keV and to have a mass ratio equal to the cosmic density ratio 
(2.68 in this case). The central degeneracy factor is \ = —0.5 
implying that the central neutrino density is about 0.4 of the 
maximum density imposed by the initial phase space constraint 
(eq. 11). The accelerations near the core radius (~ 200) kpc are 
comparable to ao 



large positive values of \ would correspond to complete de- 
generacy (P u oc p v 5 ^ 3 ) and large negative values, to complete 
non-degeneracy (P„ oc p„cr„ 2 ). 

The density and the velocity dispersion formally de- 
termine the chemical potential, but there is no inde- 
pendent method, short of numerical calculations of col- 
lapse of the two- fluid mixture to estimate \- This es- 
sentially depends upon how close the actual central den- 
sity of the neutrino fluid is to the maximum permit- 
ted by the phase space constraints and the effectiveness 
of violent relaxation in equalising the velocity dispersions 
of the two fluids. There are various arguments suggest- 
ing that the true density should range between 10% and 
100% pmax dMadsen fc Epstein 1984] IKofman et al. 19961 
Nakaj ima 8c Morikawa 2007[ ) ; here I take p v ~ pmax be- 
cause this is in rough agreement with the observational re- 
sults shown in Fig. 2 (i.e., the central neutrino density is 
50% that corresponding to quantum mechanical degener- 
acy). This corresponds to x = —0.5, and in the calcula- 
tions described below, I have fixed the central value of x 
at this value. Realistically, due to random initial conditions 
for cluster formation, one would expect a dispersion in this 
parameter. 

Assuming a comparable mass (m„ = 1.9 eV) for all 
three neutrino types, a constant temperature for both fluids 
(with the velocity dispersion of neutrinos equal to that of 
the baryons), and a mass ratio of the two components equal 
to the cosmic density ratio (2.68 in this case), the structure 



Figure 4. A log-log plot of the luminosity-temperature relation 
for self consistent MOND-neutrino- baryon fluid spheres compared 
to the observations of Ikebe et al. (2002). The neutrino mass 
is taken to be 1.9 eV and the degeneracy factor is -0.5. With 
the assumptions of isothermal state, equal velocity dispersions, 
and cosmic ratio of baryons-to neutrinos the structure and X-ray 
luminosity of these objects is completely specified. 



of the object may be determined by numerical integration of 
eqs. 18 supplemented by eqs. 19-23. With these assumptions 
and constraints, the structure completely determined when 
the gas temperature is specified. Fig. 3 shows the density 
distribution for the two fluids in the case where Tk e v = 6. 
The density of the two components effectively track each 
other in the inner regions. The distributions can be described 
as a roughly constant density core extending to about 300 
kpc, followed by a rapid decrease (asymptotically p oc r~ 3 ' 5 
for both components. 

The X-ray luminosity resulting from free-free emission 
of the hot ionised gas may also be calculated for any such ob- 
ject. The resulting luminosity-temperature relation is shown 
in Fig. 4 compared to the sample of Ikebe et al. (2002) . This 
is similar to the figure shown in Sanders (2003) but there the 
dark component was added arbitrarily, with no underlying 
physics. Here the X-ray emission is that from self-consistent 
neutrino-baryon fluid spheres with only the neutrino mass 
and degeneracy factor specified arbitrarily. The results gen- 
erally agree with the observed relation; in particular the 
steepening of the relation for low temperature clusters is 
evident. 



6 NON-ISOTHERMAL (3 MODELS 

The gas density distribution in real clusters does not re- 
semble that shown in Fig. 3. The density in these isother- 
mal models falls too rapidly beyond the core to be consis- 
tent with the observed distribution of X-ray surface bright- 
ness in clusters; i.e., the power law implied by the fitted 
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Figure 5. The predicted projected temperature distribution in 
A2256 compared to that observed (de Grandi & Molendi). The 
assumed neutrino mass is 1.9 eV. The /3 model parameters are the 
mean emission weighted gas temperature is 7.2 keV. The neutri- 
nos are assumed to have a constant velocity dispersion appro- 
priate to this temperature. The solid curve shows the prediction 
when the central gas temperature is take to be 7.5 keV and the 
dashed curves show the distributions when the central tempera- 
ture is taken to be 8 keV (upper curve) or 7 keV (lower curve). 



(3 models (eq. 6) is more like r~ 2 rather than r -3 ' 5 . How- 
ever, it is now known that the gas in clusters cannot be 
described as isothermal. There appears to be a characteris- 
tic temperature profile in clusters: in non-cooling flow clus- 
ters the temperature declines by 20 or 30% between 100 kpc 
and 1 Mpc; in cooling flow clusters there is first a rapid 
rise in the inner core region (by as much as a factor of 
two) followed by the decline observed in non-cooling flow 
clusters (|De Grandi fe Molendi "2002] IPiffaretti et al. 20051 
IPratt et al. 20061 IViFhlinin et al. 2005) 

Therefore, here I take a different approach. I assume 
the density distribution of the gas is given by a /3 model 
and I solve eq. 19 to determine the temperature distribu- 
tion for the gas. The neutrino fluid is again assumed to be 
isothermal with a velocity dispersion equal to that implied 
by mean emission weighted temperature of the gas. In this 
way, one may predict the temperature distribution in a par- 
ticular cluster with a fitted /3 model. The only adjustable 
parameter is the central value of the gas temperature, and 
this is chosen to match either the mean emission-weighted 
temperature of the gas or to achieve the best fit to the pro- 
jected gas temperature as a function of radius. In fact, there 
is generally a narrow range of this parameter which yields 
a reasonable temperature distribution for the cluster- one 
which does not increase rapidly to absurdly high values or 
one which does not fall rapidly to zero. 

This is illustrated in Fig. 5 which is the radial depen- 
dence of projected gas temperature for the model fit to 
A2256 with m v = 1.9 eV and x = —0.5. This is a non-cooling 



flow cluster with a low central density (n e = .0031 cm - 3) 
and a large core radius (r c = 419 kpc) . The decline of surface 
brightness implies (3 = 0.914 QReiprich fc Bohringer 2002p . 
Taking a central gas temperature of 7.5 keV produces a 
projected temperature profile (solid curve) which is consis- 
tent with that observed (|De Grandi fc Molendi 2002p . The 
dashed curves illustrate the effect of taking a central gas 
temperature one-half keV higher or lower; it is evident that 
the resulting temperature profile is strongly dependent upon 
the initial assumed central gas temperature. This procedure 
has been carried out for several non-cooling flow clusters, 
and the characteristic temperature profile is as shown here; 
i.e., it is consistent with the observed profiles for non-cooling 
flow clusters. 

For cooling flow clusters, the results are even more sen- 
sitively dependent upon the assumed central temperature, 
but, in general, the calculated temperature profiles agree 
with those observed for cooling flow clusters. This is shown 
in Fig. 6 which illustrates the predicted temperature profile 
(again with m v — 1.9 eV and x = —0.5) for the cooling flow 
cluster A85 with /3 model parameters of [3 = 0.532, r c = 58.1 
kpc, n e = 0.0204 cm" 3 ( |Reiprich fc Bohringer 2002} . This is 
compared to the observed temperature profile by De Grandi 
& Molendi. The dotted curves show the effect of increasing 
or decreasing the assumed central gas temperature by 0.1 
keV; i.e., the results here are extremely dependent upon the 
this parameter. 

For clusters with f3 model fits characterised by a small 
core radius and high central electron density, this is the gen- 
eral pattern predicted the two fluid models: a rapid rise in 
temperature followed by a gradual decline. In general, the 
predicted central temperature is lower than observed and the 
detailed agreement is less impressive than for non-cooling 
flow clusters. It should be kept in mind, however, that no 
baryonic component other than the gas is included in these 
calculations; the galaxies, and in particular, large central cD 
galaxies are not part of the mass modelling. 

Overall, it appears that when the j3 model is imposed 
upon the gas density distribution, the temperature distri- 
bution required by self-consistency is in general agreement 
with the observed temperature profiles. 



7 CAN INDIVIDUAL GALAXIES HAVE 
NEUTRINO DARK MATTER HALOS? 

It would seem possible that individual galaxies could possess 
extensive neutrino dark halos, even if the neutrino mass is as 
low as 2 eV. It is true that no primordial fluctuations on the 
scale of galaxies would survive in the cosmic neutrino fluid 
due to free streaming. However, in the context of MOND, 
the baryonic component would collapse first and act as a 
seed for subsequent neutrino infall. 

The neutrino halo could have a mass almost three times 
larger than the baryonic galaxy, but because of the phase 
space constraint, it would be very extensive. Rewriting eq. 
16 in terms of rotation velocity, V r , we would find, 



R c * = 1.3 





4 
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_200 km/s. 



for m„ = 2 eV the halo would extend to 500 kpc for a 
typical massive galaxy. The total mass of the neutrino halo 
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Figure 6. The predicted projected temperature distribution in 
A85 (solid curve) compared to that observed (de Grandi & 
Molendi). The assumed neutrino mass is 1.9 eV and the degener- 
acy factor is x = —0.5. The j3 model parameters are j3 = 0.532, 
r c = 58.1 kpc, n e = .0204 cm~3. The central gas temperature 
is set to 1.69 keV. The two dotted curves are the predicted tem- 
perature distributions when the central gas temperature is set to 
1.79 keV (upper curve) and 1.59 keV. This illustrates the extreme 
sensitivity of the predicted distribution to central gas tempera- 
ture in cooling flow clusters (high central electron density, small 
core radius). 

on galactic scale (10-20 kpc) would be less than 1% of the 
galaxy's baryonic mass and the contribution of the halo to 
the observed rotation curve in the outer regions would be 
less than 10 km/s. Therefore, such a halo could have no 
effect on observed galaxy kinematics. 

It would, however, affect the weak lensing properties 
of galaxies on large scale. On the scale of hundreds of 
kpc, the halo mass would exceed the mass of the galaxy, 
and dominate the lensing signal. If, for example, the halo 
were flattened, then, with sufficient statistics, this could 
be apparent from the pattern of background galaxy images 
flHoekstra, Yee fc Gladders 2004| >. 



8 CONCLUSIONS 

Non-baryonic dark matter in the form of active neutrinos 
certainly exists; only its contribution to the total mean cos- 
mic density is in question. If the mass of the three neutrino 
types is as large as 2 eV, then neutrinos will be a constituent 
of massive rich clusters of galaxies. For this mass scale, the 
cosmological density ratio is of the order of the remaining 
mass discrepancy in clusters calculated via MOND; i.e., such 
neutrinos would complete the mass budget of clusters in the 
context of MOND. Maintaining the cosmological mass ra- 
tio in clusters leads naturally to a dependence of central 
electron density on temperature (n e oc T 3 / 2 )- a dependence 
which is consistent with density inferred from f3 model fits 



to the X-ray intensity distribution in non-cooling flow clus- 
ters (Fig. 2). This gas density-temperature dependence cor- 
responds precisely to the constant entropy floor that has 
been proposed to account for the steepness of the observed 
luminosity-temperature relation, but in this case no preheat- 
ing of the cosmic baryon fluid is needed; the floor is pro- 
vided by the phase space constraints on neutrino density. 
Of course, in cooling flow clusters one would expect there to 
be an inward flow of gas and breaking of the cosmological 
density ratio- with resulting higher electron densities. 

Assuming that the cosmic density ratio of the two fluids 
is maintained and that the neutrinos and baryons have the 
same velocity dispersion, then only the neutrino mass and 
the degeneracy factor need be specified in order to calcu- 
late hydrostatic two fluid models of clusters. Such models, 
with m v « 2 eV exhibit the observed X-ray luminosity- 
temperature relation- including a break to a steeper rela- 
tion for lower temperature clusters (Fig. 4). If the gas den- 
sity distribution is constrained to follow the fitted models 
for clusters, then the temperature distribution of the gas re- 
quired for self-consistency resembles that in actual clusters- 
a general decline over hundreds of kpc for non-cooling clus- 
ters, but a rapid rise followed by a decline for cooling flow 
clusters (Figs. 5 & 6). 

Thus these calculations support the suggestion that the 
missing component in clusters of galaxies, in the context 
of MOND, may be neutrinos with mass near the present 
experimental upper limit. A remaining observational prob- 
lem with this is that some clusters- primarily cooling flow 
clusters- apparently require a higher central density of dark 
material than is permitted by the phase space constraints on 
neutrinos (Sanders 2003). This is evident both from X-ray 
observations and observations of strong gravitational lens- 
ing in the centres of some clusters where the implied central 
mass within an Einstein ring radius (100-200 kpc) may be in 
excess of 10 13 M (|Sanders 19990 . Of course, the presence 
of a massive central galaxy and its effect on the distribution 
of neutrinos has not been included here, but it may well be 
that an additional undetected baryonic component in cool- 
ing flow clusters- perhaps resulting from the cooling flow 
itself- is required. I have also not considered the possibility 
that there is at least one massive sterile neutrino- a pos- 
sibility with some theoretical and experimental motivation 
( [Bilenky et al. 20030 . In any case, the greatest part of the 
discrepancy in clusters would be accounted for by the three 
types of active neutrinos if their mass is near 2 eV. 

We will not have to wait long for this possibility to be 
falsified (or confirmed). Currently planned decay experi- 
ments ( [Osipowicz et al. 20010 will push the upper limit on 
the electron neutrino mass to a few tenths of an electron 
volt within a few years. If it turns out to be the case that 
m v m 2 eV, then, with MOND, the old problem of clusters 
is solved. 
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